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. _■ _ _ INTRODUCTION _ _ _  _ _ 

As  semiconductor  technology  continues  to  pursue  the  scaling-down  of 
IC  device  dimensions  into  the  submicron  (less  than  ten  thousand  Angstroms) 
regime,  many  novel  and  interesting  questions  will  emerge  concerning  the 
physics  of  charged  particles  in  semiconductors.  One  of  the  more  important 
topics  to  be  considered  is  that  of  the  appropriate  transport  (1)  picture 
to  be  used  for  a  given  spatial  and  temporal  regime.  Moreover,  from  the 
point  of  view  of  device  physics,  it  is  most  desirable  to  have  a  microscopic 
description  of  semiconductor  transport  which  is  computationally  manageable 
or  at  least  amenable  to  phenomenological  treatment  so  that  its  properties 
can  be  meaningfully  incorporated  into  device  simulations.  In  this  paper 
attention  is  focused  on  a  useful  transport  methodology  for  the  ultra¬ 
submicron  regime.  At  present,  military  electronic  devices  with  transit  ; 
lengths  in  the  ultra-submicron  region  are  coming  to  fruition  due  to  the 
advent  of  MBE  (2)  processing  so  that  the  need  for  an  appropriate  descrip¬ 
tion  of  transport  in  this  regime  is  indeed  imperative. 

Semiconductor  transport  in  the  ultra-submicron  regime  approaches  the 
category  of  quantum  transport.  This  is  suggested  by  the  fact  that  within 
the  effective  mass  approximation  the  thermal  deBroglie  wavelength  for 
electrons  in  semiconductors  (see  Fig.  1)  is  of  the  order  of  ultra¬ 
submicron  dimensions.  Whereas  classical  transport  physics  is  based  on  the 
concept  of  a  probability  distribution  function  which  is  defined  over  the 
phase  space  of  the  system,  in  the  quantum  formulation  of  transport  physics, 
the  concept  of  a  phase  space  distribution  function  is  not  possible 
inasmuch  as  the  non-commutation  of  the  position  and  momentum  operators 
(the  Heisenberg  uncertainty  principle)  precludes  the  precise  specification 
of  a  point  in  phase  space.  However,  within  the  matrix  formulation  of 
quantum  mechanics,  it  is  possible  to  construct  a  "probability"  density  ; 
matrix  which  is  often  interpreted  as  the  analog  of  the  classical  distri¬ 
bution  function. 
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Figure  1.  Thermal  deBroglie 
wavelength  vs.  effective  mass 


There  is  yet  another  approach  to  the  formulation  oif  quantum  transport 
based  on  the  concept  of  the"  Wigner  distribution  function  (WDF)  (3).  This 
formalism  is  particularly  attractive  for  use  in  ultra-submicron  device 
transport  in  that  it  contains  all  of  the  quantum  mechanical  information 
about  the  state  of  the  system  yet  has  elements  of  the  classical  picture 
implicitly  built  in.  Thus,  the  intent  of  this  study  is  to  explore  the 
potential  usefulness  of  the  WDF  as  well  as  other  possible  quantum  distri¬ 
bution  functions  for  describing  quantum  (ultra— submicron)  device  transport 
in  semiconductors.  To  this  end  we  first  review  the  salient  features  of 
the  WDF  and  then  discuss  a  new  result,  the  derivation  of  the  first  three 
quantum  moment  equations  using  the  WDF.  It  is  shown  that  the  moment  equa¬ 
tions  contain  quantum  corrections  to  the  classical  moment  equations;  these 
quantum  terms  are  non-negligib.le  when  the  transit  lengths  are  of  the  order 
of  the  carrier  deBroglie  wavelength. 

THE  WIGNER  DISTRIBUTION  FUNCTION 

The  Wigner  distribution  function  (3)  is  generally  defined  in  terms  of 
all  the  generalized  coordinates  and  momenta  of  the  system  as 

PW(xl*-*V  pl***-pn>  =  '^'n  dyl***dyn  T  (X1  +  T»”*Xn 
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However,  for  simplicity  in  this  paper,  we  discuss  the  properties  of  a _ 

single  coordinate  and  momentum  WDF: 

.*  \ 

PW<X*P)  “  2ih  f  dy  y*  (x  +  f)y  (x  -  |)  elpy/TS  (2)  ! 

( 

where  y(x)  represents  Jt.he_sjtate  of  the  system  in  the  coordinate  represents-, 
tion.  Although  we  will  be  treating  the  WDF  for  the  special  case  of  pure  ) 
states,  the  adaptation  of  this  formalism  to  include  mixed  states  is  accom-  j 
plished  through  the  generalization  _  ! 


Pw<x,p) 


z  p  r — i— 

n  nl™ 


J  dy  y*  (x  +  %)  Vn  (x  -  J)  eipy/1* 


where  P_  is  the  probability  to  be  in  state  "n".  (For  example,  for  a  system 

"  ~*-Eo7kY  * 

in  -contact --wi-t-h-a-  heat  -bafch-at— eonstant  temperature,  P  =  e  '  .) 

The  distribution  function  of  Eq.  (2)  has  interesting  properties  in  j 
that  the  integration  of  this  function  over  all  momenta  leads  to  the  proba-  : 
bility  density  in  real  space;  conversely,  the  integration  of  this  function  j 
over  all  coordinates  leads  to  the  probability  density  in  momentum  space.  , 
In  mathematical  terms,  ; 

|  00  1 

!  /  Pw(x,p)  dp  =  y  (x)y(x)  (3a)  I 


J  Pw(x,p)  dx  =  $*(p)<Kp), 


where 


(2irtl)~1/2  /  y(x)dx. 


It  follows  limnediately  from  Eq.  (3)  that,  for  an  observable  W(x,p)  j 
which  is  either  a  function  of  momentum  operator  alone  or  of  position 
operator  alone,  or  any  additive  combination  therein,  the  expectation  value  1 
of  the  observable  is  given  by  j 


(w)-  //  WPw(x,p)  dx  dp. 


which  is  analogous  to  the  classical  expression  for  the  average  value.  ! 

Herein  lies  the  Interesting  aspect  of  the  Wigner  distribution  function;  the 
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result  of  Eq.  (4)  suggests  that  it  is  possible  to  transfer'  many  of  the  re- 
suits  of  classical  transport  theory  into  quantum  transport  theory  by  simply 
replacing’ the  classical  distribution  function  by  the  Wigner  distribution 
function.  However,  unlike  the  density  matrix,  the  Wigner  distribution  I 
function  itself  cannot  be  viewed  as  the  quantum  analog  of  the  classical  I 
distribution  function  since  it  is  generally  not  positive  definite  and  non-  j 


unique  [Pjj(x,p)  of  Eq.  (2)  is  not  the  only  bilinear  expression  (3-5)  in  T  j 
that  satisfies  Eq.  (3)].  ! 

Further  resemblance  of  the  Wigner  distribution  function  to  the  classi¬ 
cal  distribution  function  is  apparent  by  examining  the  equation  of  time  ; 
evolution  for  Py(x,p).  Upon  assuming  that¥(x)  in  Eq.  (2)  satisfies  the  j 
Schrodinger  equation  for"  a  system  with  Hamiltonian  H  =  p%/2m  +  V(x),  it  can- 
be  readily  shown  that  Pjj(x,p)  satisfies  “the  equation  j 


3PwA»t-+  (  p/m>  SP^Ax— +-  0  ---P^  ~~  0,- 


— (5)_ 


where 


0-Pw  =  - 


2  "  .  .n  (h/2)2n+1  *2n+1' 

'Jo  <2"+X>!  3X2" 


+  1) !  3x2„+1  Sp2„+1 


It  is  evident  that  in  the  limit  h  -*■  0,  0  •  Py  in  Eq.  (6)  becomes 


(6) 


0  -Pw  -  -OV/3x)OPw/3p) 


(7) 


so  that  Eq.  (5)  reduces  to  the  classical  collisionless  Boltzman  equation.  j 
The  Wigner  distribution  function  defined  in  Eq.  (2)  is  derivable  (8)  j 
from  the  Fourier  inversion  of  the  expectation  value  (with  respect  to  state  j 

-  T(x))of  the  operator  ei(rP+6x^ (here„[x  and  p]satisfy  the  commutation  j 

relation  [x,p]  -  1ft).  As  such,  j 

j 

Pw(x,p)  ’‘-—H  Cw(T,e)e“i(Tp+0x)dT  d8,  (8a)  j 

4  IT  I 

where 

C^t.e)  -  /  T*(x)e1(rP+0x)  »  (x)  dx,  (8b) 

and  the  interval  of  integration  is  (-*,*)  unless  otherwise  specified.  In 
! order  to  show  that  the  right-hand  side  of  Eq.  (8a)  is  indeed  the  Wigner 
distribution  function  as  defined  in  Eq.  (2),  note,  from  the  Baker-Hauadorff; 
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theorem  (9) y  that  can  be  written 


ei(xp+9x)  _  eirp/2ei6xeiTp/2 


in  which  case  of  Eq.  (8b)  becomes 

CO 

C^T.e)  =  /  [e 


-itp/2  ,  (x)]  *  ei6x  [e-ixp/2  f  (x)]  dx 


(9) 


(10) 


which  further  reduces  to 

oo 

Cw(t,6)  =  /  S'  (x  -  *sxh)ei^X'i’  (x  +  *jxh)  dx..  (11) 


Then,  by  inserting  C^TjG)  of  Eq.  (11)  into  the  right-hand  side  of  Eq.  (8a),- 
integrating  over  the  variable  9  by  using  the  relation 

i  i 

7  ei0(x,“x"  >  d0  -  2ird(x’-x”),  j 

; 

> 

and  letting  x  =  y/Tt,  the  desired  result  is  obtained.  j 

The  method  outlined  above  to  arrive  at  the  Wigner  distribution  func¬ 
tion  is  based  on  the  notion  of  a  characteristic  function.  The  character-  \ 
istlc  function  of  an  observable.  A,  with  respect  to  state  (here,  the 
Dirac  notation  is  utilized  for  purposes  of  generality)  is  defined  as 

i 

A 

CA(0  -  <'T|eUA|'F),  (12)  1 

j 

where  5  is  a  real  parameter.  Assuming  A  to  possess  an  eigenvalue  spectrum  i 
given  by  a|a^  -  A*  (A’)  ,  CA(0  can  be  evaluated  in  the  A' -represen tat ion  as  j 

,  i 

i 

CA(0  -  /  dA'  /  dA"  < y | A*>  <A,|eiU|A,,>  <A"|<F>  .  (13)  j 

( 

«  f  j 

Since  (a'  |ei^A|A")  »  e*^  6(A'  -  A”)  in  the  A' -representation,  CA(£)  in  j 

Eq.  (13)  reduces  to  j 

CA(C)  “  /  dA’e15A'|TA»|2,  <*4>  J 

!  i 

where  |f^t|  •  KA’|f)J  =  P(Af),  the  probability  distribution  function  for! 
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measuring  A'  while  in  state  |*?^.  Hence,  the  characteristic  function  for 
A  is  the  fourier  transform  of  the  probability  distribution  function  P(A’).  ■ 
Subsequent  inversion  of  Eq.  (14)  above  leads  to 

{ 

P(A’)  =  /  CA(Oe“iCA'  d£..  (15)  : 

i 

The  Wigner  distribution  function  was  derived  by  taking  the  Fourier 

transform  of  the  characteristic  function  for  e*(TP+®x).  in  vfew  Qf  the 
connection  between  the  probability  distribution  function  and  the  character-! 
istic  function  for  a  given  observable,  this  approach  seems  to  be  a  natural 
way  of  obtaining  a  distribution  function  for  momentum  and  position.  Unfor-; 
tunately,  the  noncommutative  nature  of  the  two  observables  destroys  the  ; 
convenient  probability  interpretation  of  the  characteristic  function 

implicit  in  Eq.  (15)  . _ _ _ _ _ j 

i  In  order  to  demonstrate  this  point, ^assume  the  characteristic  function] 
for  two  noncommuting  observables,  A  and  B,  to  be  j 


CAbM2>  =  <'F|e±(SlA+C2B>  jy> 


Observables  A  and  B  are  assumed  to  have  eigenvalue  spectra 

i 

a|a')  -  A’ |a’}  ,  b|b’>  -  B* |b’>  . 

and  are  chosen  so  that  [A,  [A,B]]  =  [B,  [A,B]]  =  0. 

This  assumption  is  imposed  so  that  the  identity 

ei(C1A+C2B)  „  eiexA  eiC2B  e~5lMA>®J/2 


(17)  i 


(18)  | 


may  be  used.  j 

Inserting  Eq.  (18)  into  Eq.  (16)  while  obtaining  the  matrix  elements  j 

of  ei^lA  in  the  A-representation  and  e1^®  in  the  B-representation  results  : 
in 

CAB<CiC2)  -  e“5lMA»®  1/2  /  dA'  /  dB’  ei( 5 1  A*+e2®  ’ (A'Xa'  |b'><B*  |v>(19)  \ 


>In  Eq.  (19),  it  is  assumed  that  [A,B]  is  a  c-number  independent  of  the 
eigenvalues  A'  and  B'.  We  define  E^A'.B’),  the  generalized  Wigner  distri¬ 
bution  function,  to  be 
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'  F(A*  ,B')  =  .<» I A*XV  , 


(20) 


so  that  “ 


F(A'.B’) - l-^r  /dC  JdE  eClC2fA‘Bl/2  CA_(ei#e2)e“1<€lA,+52B,).  (21) 

(2iOZ  12  AB  i  ^ 


—ft —is- -evident  from  Eqs.  (20,21)  that 


/  F(A’.B')  dA*  =-|<B'|¥>}2  =  ^  f  d52  ^^(O.^e-1^8' 


(22a) 


and 


J  F (A* , B * )~dB ’  =  KaM*>|2  -  “  /"d?!  yS^O)#'15/'.  (22b)  ! 


Thus,  Eq.  (21)  establishes  the  relationship  between  the  characteristic 
function  for  two  arbitrary  noncomauting  observables  and  the  generalized  ! 
Wigner  distribution  function.  The  generalized  distribution  function  has  , 
the  essential  properties  of  the  conventional  Uigner  function  in  that  an 
integration  of  the  generalized  function  over  the  eigenvalue  spectrum  of  one 
observable  leads  to  the  probability  density  in  the  canonically  conjugate 
observable  [  Eq.  (22)].  j 

There  is  no  simple  probability  interpretation  of  F(A',B*)  in  Eqs.  (20* 
21)  because  of  the  necessary  overlap  between  the  states  of  the  noncommuting 
observables.  However,  if  A  and  B  are  made  to  commute  so  that  |A'^  and  |b')  | 
are  a  common  set  of  eigenvectors,  then  F(A',B*)  reduces  to  the  probability: 
distribution  function  for  A  and  8. 


Finally,  it  is  ngted  that  the  conventional  Wigner  distribution  func-  > 
tion  for  observables  A  and  B  is  obtained  from 


'71^  UC‘  ^  C*B<51^2)e"1<5lA’+{2B  )- 


(23) ; 


■with  C^b(?iC 2)  defined  in  Eq.  (16),  whereas  the  alternative  distribution  1 
function,  F(A',B’),  introduced  in  Eqs.  (20,21)  differs  from  the  Wigner 

function  due  to  the  presence  of  the  phase  factor  e^1*2  *  in  the  inte-  ’ 
grand  of  Eq.  (21).  For  A  ■  x  and  B  *  p,  fy(x,p)  in  Eq.  (23)  reduces  to  the 
Wigner  function  of  Eq.  (2),  whereas  F(x,p)  defined  from  Eq.  (20)  becomes 


CO 

F(x,p)  -  J  dy  Y*  (x)'F(x  -  y)eipy/1i  =(2nh)"1/2  (x)elpx/Tk  *(p),  (24) 


IAFRATE 


where  <|>(p)  'is  defined  in  Eq.  (3b).  It  is  evident  that  there  is  a  family 
of  functions  which  are  bilinear  in  ’?  yet  satisfy  the  sum  rules  of  . 

Eqs.  (3a ,b). 

There  are  some  interesting  questions  to  be  resolved  concerning  the 
uniqueness  and  positive  definiteness  of  Wigner-type  quantum  distribution 
functions.  Nevertheless,  these  distribution  functions  serve  a  useful  pur- 
pose  for  calculating  quantum  mechanical  observables  in  transport  (7) 
studies  and  numerous  solid-state  (10,11)  problems ~ 

■-  MOMENT  EQUATIONS 

In  this  section,  we  derive  the  firsF  "three  moments  of  a  "Wigner- 
Boltzmann"-like  transport  equation. 


aV  p  jp* 

at. m-  -ax— ~  -  -V _ I.  at. 


y 


coll. 


-  ..~(2s)~i 


'This  equation  was  constructed  to  include  an  ad-hoc  collision  term  which  may 
not  necessarily  express  the  same  phenomenology  as  that  of  the  classical 
Boltzmann  transport  equation,  since  Pw  is  not  a  true  probability  distritu-  ; 
tion  function.  These  problems  are  conceptually  reduced  when  dealing  with 
moments  in  a  relaxation  approximation.  The  moment  equations  are  obtained  f 
by  multiplying  Eq.  (25)  by  an  appropriate  function  of  momentum,  4>(p),  and 
then  integrating  over  all  momenta  to  obtain: 


3<$>  ,t.  1  3 

3t  m  3x 


«|>p> 


-  E 
n=0 


(~) 

2i 


Tttiji 


.  2n+l 


3x 


2n+l 


32n+l 

v(x)  /  Hp)  P 


(26) 


coll. 


where  <  >  refers  to  an  integration  over  momentum.  In  making  the  assump¬ 

tions  that  <Kp)  be  an  analytic  function  of  momentum  and  that  Pw(x,p)  vanish 
at  the  momentum  limits,  it  follows  that 


/  <Kp)  ( 


d2n+1 

3p2n+1 


V  ^ 


<-l)2“+1  / 


3p2"+l 


<Kp)>  Pwdp 


(27) 
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and  that  Eq.  (26)  becomes 

+ 1  ( A“  -fa  (£=L  « <  i=2f 

n=0  2i  <2n+1>J. . 3x2n+l .  . . 3p2n+l. 


na<  ^  >  -  <  ♦  O  > 

. ap2n+!  -  at  con. 


'Z 

For  specific  values  of  4>(p)  =  p  ,  P»..  a?4  .‘2m»  E9*  becomes 

£  +  i  A<^-«, 


(29a) 


- <29b)! 


-L  a  <p2>+  1  S  <3>+.<Ii  M.<£i  ("£)  >  (29c) 

2m  3t  v  2mz  3x  m  3x  2m  9t  . . 


coll. 


where  upon  reduction 


<Pn>  -  (“)  £  (-D 


j  nt  9JV  (a 
(n~j)!j!  3xj 


a  Jy(x) 

}  * 


thereby  showing  the  dependence  of  <pn>.  For  n  =  0,  1,  2,  3  we  show  <p  > 
explicitly! 


<P°>  =  f*» ,  <P1>-^(^S-»^T). 


Ola) 


<p2>  .  ( JL)2  (y*  if*  -  2  111  21  +  l£l%) 

<P  >  to7  3x  3x  3x  '» 


(31b)  ; 


^  (V*  a3y  7  34'*  ill  +  3  w  _  ili  V)  nic) 

<P3>  -  (77)  (*  a^-3-sr  3^+3~3^  33?  Ulc) 


Note  that  use  of  Eqs.  (31a)  in  Eq.  (29a)  results  in  the  correct  quantum 
mechanical  continuity  equation. . .  _ _ _ _ _ 
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In  order  to  see  the  transition  between  the  quantum  and  classical 

regimes,  ve  invoke  a  wavef unction  of  the  form  ¥(x,t)  =  A(x,t)e^^,t:^^, 
and  so 

A2(x,t)  =  ¥*f  =  p(x,t)  ,  ~  3-S3(xX>t)  =  v(x,t)  (3 


where  p(x,t)  is  the  probability  density  and  v(x,t)  is  the  ensemble  veloc 
ity.  Inserting  the  wavefunction  of  Eq.  (32)  into  Eq.  (31)  results  in 


It  is  evident  from  Eq.  (33)  that  the  terms  possessing  an  explicit 
dependence  on  "ft"  are  the  quantum  corrections  to  the  classical  momentum- 
density  moments.  Furthermore,  insertion  of  these  momentum-moments  into 
Eqs.  (29)  results  in  a  set  of  moment  equations  which  contain  explicit 
quantum  corrections  as  well.  Of  course,  as  lT-»-0,  these  moment  equations 
reduce  to  the  classical,  zero  temperature,  moment  equations. 

The  logarithmic  derivative  term  appearing  in  Eqs.  (33)  is  not  uncom¬ 
mon;  such  a  term  appears  in  the  real  part  of  the  Schrodinger  equation  when 
the  wavefunction  associated  with  Eq.  (32)  is  used  (9).  In  the  limit  where 
this  term  is  negligible,  the  Schrodinger  equation  goes  over  to  the 
Hamilton-Jacobi  equation. 

We  have  made  order  of  magnitude  estimates  of  the  strengths  of  the 
quantum  correction  term  appearing  in  Eq.  (33).  Using  a  Gaussian  spatial 
variation  for  p(x),  and  a  thermal  ensemble  value  for  mv  ,  typical  of 
central  valley  GaAs  electrons,  we  find  that  the  quantum  correction  is  sub¬ 
stantial  for  distances  of  the  order  of  100  A. 


\  SUMMARY 

In  this  paper^weT hav<^ puts  forth  a  formalism  for  treating  ultra- 
submicrofi  device  transport.  The  formalism  results  in  a  useful  and 
attractive  methodology  for  describing  ^quantum*- device  transport  in  that 
the  theory  is  derived  from  a  fully  quantum  mechanical  representation,  yet 
implicitly  contains  elements  of  the  semiclassical  semiconductor  transport 
picture.  .  The  basic. three  semiconductor  quantum  transport  equations  were 
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"derived  using  the  Wigner  distribution  function,  these  transport  equations 
were  shown  to  contain  explicit  quantum  corrections;  these \quahtum  correc¬ 
tions  are  non-negligible  when  the  transit  lengths  of  the  semiconductor 
device  are  of  the  order  of  the  carrier  deBroglie  wavelength.  Since  the 
carrier  deBroglie  wavelength  for  carriers  is  of  the  order  of  hundreds  of 
angstroms  in  III-V  device  materials  of  interest  to  the  Army,  the  quantum 
description  of  transport  as  described  here  will  play  a  vital  r0le  in 
predicting  the  electrical  behavior  of  present  and  future  generation  ultra- 
submicron  devices.  \ 

In  future  studies,  wa  plan  to -expand -ear  theoretical  effort&^to 
include  device  modeling  of  submicron  and  ultra-submicron  semiconductor 
devices  such  as  P-N  junctions,  planar  doped  barriers,  and  one-  and  two- 
dimensional  superlattices;  for  these  devices,  it  is  clear  that  quantum 
transport  will  indeed  be  necessary  to  explain  their  semiconductor  transport 
characteristics.  _  __  _  _ 
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